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This paper reports results for directed flow vi and elliptic flow V2 of charged particles in Cu+Cu 
collisions at ^/snn = 22.4 GeV at the Relativistic Heavy Ion Collider. The measurements are for 
the — 60% most central collisions, using charged particles observed in the STAR detector. Our 
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measurements extend to 22.4 GeV Cu+Cu collisions the prior observation that vi is independent 
of the system size at 62.4 and 200 GeV, and also extend the scaling of Vi with r;/?/bcam to this 
system. The measured V2(pr) in Cu+Cu collisions is similar for ^/snn throughout the range 22.4 
to 200 GeV. We also report a comparison with results from transport model (UrQMD and AMPT) 
calculations. The model results do not agree quantitatively with the measured vi{rf), V2{pt) and 
v 2 (v)- 

PACS numbers: 25.75.Ld, 25.75.Dw 



I. INTRODUCTION 

The study of collective flow in relativistic nu- 
clear collisions has potential to offer insights into 
the equation of state of the produced matter [l|, Q . 
Anisotropic flow is conveniently characterized by the 
Fourier coefficients Q 

v n = (cosn(0- (f) 

where the angle brackets indicate an average over all 
the particles used, <f> denotes the azimuthal angle of 
the outgoing particles, is the orientation of the 
reaction plane, and n denotes the harmonic. The 
reaction plane is defined by the beam axis and the 
vector connecting the centers of the two colliding 
nuclei. The estimated reaction plane is called the 
event plane, and its orientation is denoted above by 
^r. The procedure used in the present study to 
estimate this angle is explained in Section II. C. 

Directed flow, v±, is the first harmonic coeffi- 
cient of the Fourier expansion of the final-state 
momentum-space azimuthal anisotropy, and it re- 
flects the collective sidewards motion of the particles 
in the final state. Both hydrodynamic and nuclear 
transport models [1, [f| indicate that directed flow 
is a promising observable for investigating a possi- 
ble phase transition, especially in the region of rel- 
atively low beam energy under investigation in the 
present paper @, 0j. In particular, the shape of v\ 
as a function of rapidity, y, in the midrapidity re- 
gion is of interest because it has been argued that it 
offers sensitivity to crucial details of the expansion 
of the participant matter during the early stages of 
the collision. The models indicate that the evolving 
shape and orientation of the participant zone and 
its surface play a role in determining the azimuthal 
anisotropy measured among these particles in the 
final state. 

It has been known for a long time that in the gen- 
eral vicinity of beam rapidity, the directed flow is 
dominated by a "bounce-off" effect, whereby par- 
ticles are preferentially emitted within the reaction 
plane, and are collectively deflected towards the side 
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that corresponds to being repelled from the partic- 
ipant zone Hi]. At RHIC energies, where the ra- 
pidity gap between beams is large (e.g., more than 
ten units at ^/snn = 200 GeV), the mid-rapidity 
region cannot be treated in terms of a monotonic in- 
terpolation between the rapidity regions on either 
side. Over a region of ^/snn spanning ~ 20 to 
200 GeV, it is inferred that the slope dv\ /d/q (where 
rj is pscudorapidity) exhibits the algebraic sign as- 
sociated with bounce-off near beam rapidities, but 
then v\ (rj) crosses zero at points well away from mid- 
rapidity, and the slope dvi/drj near mid-rapidity has 
the opposite sign [Io|, [Tl| ■ 

The above phenomenon of dvi/drj or dv\/dy hav- 
ing opposite sign near mid-rapidity had been pre- 
dicted by hydrodynamic and nuclear transport mod- 
els ^2Hla | and variations on this theme had been 
given names like "anti-flow" , or "third flow compo- 
nent" or "wiggle". It has been argued that this is 
a possible signature of a phase transition between 
hadronic matter and quark gluon plasma, especially 
if it is observed for identified baryons |15j . How- 
ever, it is also possible to explain the qualitative 
features of the sign reversal in dv\/dy in a purely 
hadronic picture by assuming strong but incom- 
plete baryon stopping, together with strong space- 
momentum correlations caused by transverse radial 
expansion [l4j]. 

Elliptic flow, U2, is the second harmonic coefficient 
of the Fourier expansion. The initial-state spatial ec- 
centricity of the participant zone drives the process 
whereby the interactions produce an anisotropic dis- 
tribution of momenta relative to the reaction plane. 
The elliptic anisotropy saturates quite early in the 
collision evolution, although a little later than when 
directed flow saturates [all- Elliptic flow can pro- 
vide information about the pressure gradients in a 
hydrodynamic description, and about the effective 
degrees of freedom, the extent of thermalization, 
and the equation of state of the matter created at 
early times. Studying the dependence of elliptic flow 
on system size, number of constituent quarks, and 
transverse momentum or transverse mass, are cru- 
cial to the understanding of the properties of the 
produced matter [l[. 

We find that directed flow violates the "entropy- 
driven" multiplicity scaling which dominates all 
other soft observables. STAR has reported an in- 
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triguing new universal scaling of the phenomenon 
with collision ccntrality [ll[ . Unlike the ratio of the 
elliptic flow parameter to the system initial ec- 
centricity (e), which scales with the particle density 
in the transverse plane, v\(r\) at a given central- 
ity is found to be independent of the system size, 
and varies only with the incident energy. The dif- 
ferent scalings for Vile and v\ might arise from the 
way in which they arc developed: to produce V2, 
many momentum exchanges among particles must 
occur, while to produce V\, an important feature of 
the collision process is that different rapidity losses 
need to occur for particles at different distances from 
the center of the participant zone. This later quan- 
tity is related to incident energy. A recent (3+1)- 
dimensional hydrodynamic calculation has success- 
fully reproduced various aspects of the measured di- 
rected flow for Au+Au and Cu+Cu at 200 GeV [3. 
It is of interest to see if these features are followed at 
lower beam energies, as investigated in the current 
work. 

We report here the first measurements of the di- 
rected and elliptic flow in Cu+Cu collisions from the 
STAR [l7| experiment at the Relativistic Heavy Ion 
Collider (RHIC) at = 22.4 GeV. The mea- 

surement of directed flow is presented as a function 
of pseudorapidity (77) and elliptic flow is presented as 
a function of pseudorapidity, transverse momentum 
(pr) and centrality. 

The paper is organized as follows: Section II 
briefly describes the detectors used and provides de- 
tails of the analysis methods. In Section III, we 
present results on directed and elliptic flow (v\ and 
V2) and compare these with models. A summary is 
provided in Section IV. 



II. EXPERIMENT AND ANALYSIS 
A. STAR detector subsystems 

The results presented here are based on data 
recorded by the STAR detector [l7| . The Time Pro- 
jection Chamber (TPC) [18| is the primary tracking 
device at STAR. It is 4.2 m long and 4 m in diameter 
and its acceptance covers ±1.0 units of pseudorapid- 
ity and has full azimuthal coverage. The charged 
particle momenta are measured by reconstructing 
their trajectories through the TPC. 

The charged particle reconstruction at forward ra- 
pidities is provided by STAR's Forward Time Pro- 
jection Chambers (FTPCs) 0. The two FTPCs 
are located around the beam axis, one at each end 
of the STAR detector approximately 2.3 m from the 
nominal interaction region (IR), with acceptance in 
the pseudorapidity range 2.5 < \rj\ < 4.0, and with 
full azimuthal coverage. 



The Beam Beam Counter (BBC) detector subsys- 
tem [20} consists of two detectors mounted around 
the beam pipe, each located outside the STAR mag- 
net pole-tip at opposite ends of the detector approx- 
imately 375 cm from the center of the nominal IR. 
Each BBC detector consists of nearly circular scin- 
tillator tiles arranged in four concentric rings that 
provide full azimuthal coverage. The two inner rings 
have a pseudorapidity coverage of 3.3 < \t}\ < 5.0 
and are used to reconstruct the first harmonic event 
plane (^i), for the directed flow analysis. The outer 
BBC tiles were not used. 



B. Experimental data sets 

The results presented in this paper arc from 
Cu+Cu collisions at the nucleon-nucleon center of 
mass energy of y / s]v r /v = 22.4 GeV with a minimum 
bias trigger based on BBC coincidences [2(| Hl| . We 
present the results for the 0-60% most central events, 
for which the trigger efficiency was uniform. The 
primary collision vertex position along the beam di- 
rection (V z ) has a broad Gaussian distribution with 
a root-mean-square of 62 cm. Only events within 30 
cm of the center of the detector are selected for this 
analysis. This value is chosen as a compromise be- 
tween uniform detector performance within \rj\ < 1.0 
and sufficient statistical significance of the measured 
observables. Event vertex is further required to be in 
the transverse direction within 2.0 cm from the cen- 
ter of the beam pipe in order to reject events which 
involve interactions with the beam pipe and to min- 
imize beam-gas interactions. For directed flow anal- 
ysis we further apply an event cut based on the max- 
imum value of ADC signal from the BBC to avoid 
saturation of the BBC tiles. After these event cuts 
were applied the sample used for the directed flow 
analysis contained 350k events while the sample for 
the elliptic flow analysis consisted of 800k events. 

Centrality classes in Cu+Cu collisions at ^snn = 
22.4 GeV are defined using the number of charged 
particle tracks reconstructed in the TPC within 
pseudorapidity < 0.5 and passing within 3 cm 
of interaction vertex. The uncorrected, i.e. not cor- 
rected for acceptance and reconstruction efficiency, 
multiplicity (iV c h) distribution for events with a re- 
constructed primary vertex is shown in Fig. [T] 

The directed and elliptic flow analyses were car- 
ried out on tracks that had transverse momenta 
0.1 < pr < 4.0 GeV/c, passed within 3 cm of 
the primary vertex, had at least 15 space points 
in the main TPC acceptance (\r]\ < 1.0) or 5 space 
points in the case of tracks in the FTPC acceptance 
(2.5 < |rj| < 4.0) and the ratio of the number of ac- 
tual space points to the maximum possible number 
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FIG. 1: Uncorrected multiplicity distribution with |r/| < 
0.5 in Cu+Cu collisions at ysjvjv = 22.4 GeV. Events 
with N c h > 10 are selected for the present analysis. 



of space points for that tracks trajectory was greater 
than 0.52. 



Flow methods 



At 



(5) 



The corresponding flow harmonic is here denoted as 
w n {BBC}. In all cases, the flow harmonic or coeffi- 
cient is calculated by 



(cos n(0 - vp n )) 
(cosn(#„-* R )}' 



(6) 



where is the true reaction plane angle, is 
the event plane and the angle brackets here denote 
an average over all the particles in a specific bin in 
ccntrality. Tracks used for the calculation of v n are 
excluded from the calculation of the event plane to 
remove self-correlation. The finite number of tracks 
limits the angular resolution of the reconstructed 
event plane. Consequently, the flow coefficient that 
is obtained with the reconstructed event plane, i.e. 
the numerator in equation (6), requires an event 
plane resolution correction [3[, which is the denom- 
inator in equation (6). The event plane resolution 
is estimated from the correlation of the event planes 
of two sub-events In the case of elliptic flow 



For any Fourier harmonic n, the event flow vec- 
tor (Qn) and the event plane angle (^n) are defined 

by 1 



Q„ cos(n*„) = Q nx = ^Wi cos(n0j), (2) 



Q n sin(nf r , 



Qny — ^2 Wi Sm ( n ^ 



_i ov, /rij 



(3) 



(4) 



where sums extend over all particles i used in the 
event plane calculation, and <f>i and to, are the lab- 
oratory azimuthal angle and the weight for the i-th 
particle, respectively. In the case where we construct 
the event plane using tracks reconstructed from hits 
in the TPC or FTPC, weight for each particle is 
taken to be pr in GeV/c up to 2.0 GeV/c and con- 
stant at 2.0 GeV/c for higher pt 0- I* 1 this case 
we denote the resulting flow harmonic as 1>„{TPC} 
or t;„{FTPC}. In those cases where the event plane 
is constructed using particle trajectories determined 
from hits in the BBC detectors, 4>i denotes the fixed 
azimuthal angle of the center of the ith BBC tile, 
and Wi is the energy deposition (the ADC signal, 
Ai) in the ith BBC tile and to, is calculated from 
the ADC signals Ai, where 
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FIG. 2: The event plane resolution measured using the 
TPC (second order) and using the BBC (first order) are 
shown as a function of collision centrality for Cu+Cu 
collisions at ^/snn = 22.4 GeV. Errors are statistical 
only. 

analysis, three different ways have been tested for 
choosing the sub-events: (1) Particles with pseudo- 
rapidity — 1.0 < r\ < — 0.3 are assigned to one sub- 
event and particles with 0.3 < rj < 1.0 to the other 
sub-event. The gap between the two pscudorapidity 
regions ensures that short-range correlations such as 
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Bosc-Einstcin interference or Coulomb final-state in- 
teractions contribute negligibly to the observed cor- 
relations HU, and also may reduce possible jet ef- 
fects, such as the minijet interpretation discussed in 
Ref . j23j . (2) Particles are assigned randomly to two 
sub-events. (3) Positive particles are assigned to one 
sub-event and negative particles to the other. Figure 
[5] shows the event plane resolution as a function of 
ccntrality, where the TPC second-order event plane 
resolution was determined using the first (pscudora- 
pidity) method described above. The average event 
plane resolution for the n = 2 plane for ^{TPC} is 
0.26 ± 0.01 for collisions with centrality 0-60 %. 

The BBC event plane obtained from one BBC de- 
tector is called a sub-event plane. A combination 
of the sub-event plane vectors for both BBC detec- 
tors provides the full event plane. In the «i{BBC} 
method, we used the BBC full event plane to obtain 
directed flow in the TPC pseudorapidity range (|7j| < 
1.0). A self-correlation arises if v\ is obtained us- 
ing particles from the same pseudorapidity region as 
used for the event plane reconstruction. This prob- 
lem can arise in the «i{BBC} analysis, because there 
is partial overlap in pseudorapidity coverage between 
the FTPC and the BBC. To avoid this, when Vi was 
obtained in the FTPC coverage -4.0 < ?y < -2.5, 
the event plane was constructed using the BBC cov- 
ering 3.3 < rj < 5.0, and conversely, when using the 
FTPC coverage 2.5 < r\ < 4.0 to determine V\, the 
event plane was determined with the BBC covering 
—5.0 < rj < —3.3. The first order event plane reso- 
lution as a function of centrality is shown in Fig. O 
The average event plane resolution for the n = 1 
plane for «i{BBC} is 0.16 ± 0.03 for collisions with 
centrality 0-60%. 



III. RESULTS 

A. Directed Flow Results 

Figure [3] shows charged hadron «i{BBC} in 
Cu+Cu collisions for 0-60% centrality at ^/snn = 
22 A GeV as a function of 77, compared to that for 0- 
40% central Au+Au collisions at ^/snn = 19.6 GeV 
measured by the PHOBOS experiment [T3] . The 
PHOBOS results arc quite similar, notwithstand- 
ing the difference in system size, and the fact that 
the centrality range and beam energy are not the 
same. At 200 GeV and 62.4 GeV, we have pre- 
viously reported that directed flow is not different 
within errors for Au+Au and Cu+Cu [ll|. We find 
that this behavior extends to lower energies. Di- 
rected flow provides information about the collision 
process that complements the more widely studied 
elliptic flow. Elliptic flow is developed after a num- 
ber of momentum exchanges among particles, and 
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FIG. 3: Charged hadron ui{BBC} vs. 77 for 0-60% 
centrality Cu+Cu collisions at ^snn = 22.4 GeV. The 
errors shown are statistical. Systematic errors are dis- 
cussed in Section III.C. Results are compared to vi from 
0-40% centrality Au+Au collisions at ^/snn = 19.6 GeV 
from the PHOBOS collaboration 0. 



the number of such exchanges depends on the di- 
mensions of the participant system and on its den- 
sity. Consequently, for a given collision centrality, 
elliptic flow varies with the mass of the colliding nu- 
clei. In contrast, the observation that directed flow 
does not vary with the mass of the colliding nuclei 
is a reflection of the different mechanism that gen- 
erates v\. here, the relevant feature is the rapidity 
shift undergone by particles that are initially located 
at different distances from the center of the partic- 
ipant volume 14| - a fundamental characteristic of 
the relativistic heavy-ion interaction process. 

In Fig. [4] we compare our measurements to the 
results of the A Multi Phase Transport (AMPT) 
[24l ] and Ultra Relativistic Quantum Molecular Dy- 
namics (UrQMD) [25| models. Around midrapid- 
ity, the models predict substantially smaller slope of 
V\ (77) than observed in the data, whereas at forward 
rapidities, the models differ among themselves and 
bracket the data. The fact that the tested models do 
not reproduce the observed pattern of v\ as a func- 
tion of pseudorapidity implies the need for further 
evolution in the model descriptions. 

Figure [5] shows charged hadron v\ as a function of 
pseudorapidity scaled by the respective beam rapid- 
ity (z/beam) values for the three beam energies 22.4, 
62.4 and 200 GeV in Cu+Cu collisions (ll|. For 
Au+Au collisions over a range of ^/snn spanning 
19.6 to 200 GeV, it is an empirical observation that 
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FIG. 4: Comparison of the measured ui{BBC} as a 
function of rj in 0-60% Cu+Cu collisions at y^ijvjv = 22.4 
GeV with model predictions. The inset shows the central 
r] region in more detail. The errors are statistical only. 



results reported here for Cu+Cu at 22.4 GeV extend 
the range of applicability of this scaling behavior. 
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FIG. 6: Charged hadron t?l aS cl function Of Tj - ?/beam val- 
ues for the three beam energies 22.4, 62.4 and 200 GeV. 
The results for 62.4 and 200 GeV are for 30-60% central- 
ity Cu+Cu collisions previously reported by STAR [llj |. 
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FIG. 5: Charged hadron vi as a function of r], scaled 
by the respective ybeam for the three beam energies 22.4, 
62.4 and 200 GeV. The results for 62.4 and 200 GeV 
are for 30 — 60% centrality Cu+Cu collisions previously 
reported by STAR [U. For 22.4 GeV, the plotted results 
are for — 60% centrality. 



vi(r] i 'j/beam) lies close to a single common curve for 
all beam energies [1(1 EJ, and this type of scaling 
was first observed by NA49 at the SPS [26|. The new 



Figure [S] shows charged hadron v\ as a function 
of rj - j/bcam, i.e. in the projectile frame for three 
beam energies 22.4, 62.4 and 200 GeV [ll|. The data 
support the limiting fragmentation hypothesis [ll| in 
the region -2.6 < 77 - j/boam < 0. 



B. Elliptic Flow Results 

Figure [7] shows V2(pr) for charged hadrons from 
Cu+Cu collisions at ^/snn = 22.4 GeV mea- 
sured with the subevent method with a pseudora- 
pidity gap of 0.3 units. Also shown are the pre- 
viously published STAR results for 200 and 62.4 
GeV Cu+Cu [13 measured with the full TPC event 
plane method v 2 {TPC} and full FTPC event plane 
method ^{FTPC}. The significant observed differ- 
ence between w 2 {TPC} and^FTPC} at 200 GeV 
is discussed in detail in Ref. [231 and is attributed to 
the relatively large non-flow present in V2{TPC} for 
Cu+Cu at the top RHIC energy. We observe that 
the elliptic flow at 22.4 GeV is systematically lower 
than w 2 {TPC} at 200 GeV, however, it is similar 
to w 2 {FTPC} at 200 and 62.4 GeV, consistent with 
the earlier observation [28[ . For comparison, we also 
show V2(pt) from the UrQMD and AMPT models. 
The models do not agree with the data, but they do 
show an increase in vi with transverse momentum 
similar to the data and plateau at much lower values 
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FIG. 7: Elliptic flow versus pr for charged hadrons from 
Cu+Cu collisions at 0-60% centrality at ^/snn = 22.4 
GeV measured with the subevent method with a pseudo- 
rapidity gap of 0.3 units compared with previously pub- 
lished STAR results for 200 and 62.4 GeV Cu+Cu [13 
measured with full TPC event plane method «2{TPC} 
and full FTPC event plane method u 2 {FTPC}. The er- 
ror bars are statistical. Results are also compared to 
V2(pr) model calculations. 



FIG. 8: Elliptic flow V2(rj) for charged hadrons from 
Cu+Cu collisions at 0-60% centrality at ^/snn = 22.4 
GeV. The present STAR results are compared to the 
measurement from the PHOBOS [2{| collaboration for 
Cu+Cu at 22.4 GeV. The PHOBOS results include sta- 
tistical and systematic errors whereas the STAR results 
are plotted with statistical uncertainties only, and sys- 
tematic errors are discussed in Section III C. Results are 
also compared to V2{rf) calculations from the indicated 
models. 



of pr- The small sample size in the present analy- 
sis precludes an extension of the measurements to 
identified particle v^. 

Figure |S] shows vi(r\) for charged hadrons from 
Cu+Cu collisions at 0-60% centrality at y/s^N = 
22.4 GeV. These STAR results are compared to pub- 
lished measurements from the PHOBOS Collabora- 
tion for 0-40% central collisions at ^/snn = 22.4 
GeV The PHOBOS error bars include statis- 

tical and systematic errors, whereas the STAR data 
are plotted with statistical error bars only. System- 
atic errors are discussed in Section III.C. The STAR 
results for V2(rj) are consistent within errors with the 
PHOBOS data. We also compare with correspond- 
ing predictions from the AMPT and UrQMD mod- 
els. These models underprcdict the data at midra- 
pidity, but do show a trend that is similar to the 
data for \t)\ > 2.0. 

Figure GO presents w 2 {TPC} for 0.1 < p T < 2.0 
GeV/c and \r)\ < 1.0, as a function of centrality 
for charged hadrons from Cu+Cu collisions at 0- 
60% centrality at ^/snn = 22.4 GeV, plotted along 
with previously published results from the STAR 
collaboration at 200 GeV and 62.4 GeV [27|. The 



present W2{TPC} result was obtained using r\ sub- 
events [22| with a gap of 0.3 units in pseudorapid- 
ity while the published results are based on the full 
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FIG. 9: Elliptic flow « 2 {TPC} and v 2 {4} as a function 
of centrality for charged hadrons from Cu+Cu collisions 
at — 60% centrality at ^snn = 22.4 GeV, compared 
with previously published results from the STAR collab- 
oration at 200 GeV and 62.4 GeV. 
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TPC. Also shown in Fig. [HI is a four-particle cumu- 
lant [3(J analysis 772 {4}. This analysis method helps 
reduce some types of systematic error (see below) 
but generally requires more statistics than a method 
like ^{TPC}. The statistical errors on our i>2{4} 
measurements are small enough to be useful in the 
case of the more central collisions. The results from 
this t>2{4} analysis agree, within statistical errors, 
with the f2{TPC} results, although there are hints 
of t>2{4} being systematically lower than tj 2 {TPC} 
for the more central data points. The beam energy 
dependence of the p-r-integrated tj 2 mainly comes 
from the energy dependence of the mean px and 
the difference between the event plane reconstruc- 
tion with and without a pscudorapidity gap. 



C. Systematic uncertainties 

Any analysis of collective flow needs to consider 
the possible systematic uncertainty arising from 
non-flow Q, which refers to azimuthal correlations 
not related to the reaction plane orientation. Non- 
flow can arise from resonances, jets, strings, quan- 
tum statistics effects, final state interactions (partic- 
ularly Coulomb effects), and momentum conserva- 
tion. Different methods used to measure anisotropic 
flow are affected by non-flow in different ways, and 
are used in this analysis to guide our estimates of the 
systematic uncertainty. In general, non-flow arising 
from jets can be expected to be less troublesome at 
lower beam energies. Moreover, as described earlier, 
the relatively large pscudorapidity gap between the 
STAR TPC and the BBCs is helpful in suppressing 
non-flow. 

When the event plane is determined from a de- 
tector that is not symmetric around rj = 0, we need 
to account for correlations in the measured directed 
flow due to momentum conservation 31 1. The de- 



sired BBC 77 symmetry is present for our V\ analysis 
in the 77 region of the central TPC, but is a source of 
possible concern for the FTPC 77 region. The over- 
lap in 77 acceptance between the BBC and FTPC 
is only partial, and therefore it is feasible to com- 
pare ui{BBC}f u n (where both east and west BBCs 
are used), with ui{BBC} su b, using either the east 
or west BBC event plane for 2.5 < \ii\ < 3.3. We 
find that the difference is less than 10%, and an ex- 
trapolated average correction has been applied to 
tji{BBC} for |r?| > 3.3. 

The measured v\ must be antisymmetric about 
mid-pscudorapidity within statistical errors. Any 
larger difference is due to systematic errors. Previ- 
ous detailed studies point to the maximum forward- 
backward difference as a viable estimate of the over- 
all systematic uncertainty when the pscudorapidity 
gap is large [22[. We conclude that the overall sys- 



tematic uncertainty in our determination of v\ is ap- 
proximately 15% in the FTPC region, and about 
10% in the central TPC region. 

In our analysis of elliptic flow, tj 2 , unlike for di- 
rected flow, v±, we do not have the advantage of a 
wide 77 gap to help ensure that non-flow background 
effects are minimized. To study possible systematic 
effects associated with short range non-flow correla- 
tions, a four-particle cumulant [301 ] analysis i>2{4} as 
a function of centrality has been investigated, and 
is plotted in Fig. [5] This method suppresses non- 
flow correlations involving fewer than four particles. 
The statistical errors on our tj 2 {4} measurements are 
small enough to be useful for the most central colli- 
sions, but grow to a few tens of percent at the other 
end of our studied centrality range. The u 2 {TPC} 
and i>2{4} measurements agree within statistical er- 
rors, although the observed systematic difference, 
which might arise from non-flow effects, amounts to 
about 9% for the — 10% most central collisions. 

In the elliptic flow measurement, we have used 
the 77 sub-event method with a gap of 0.3 units in 
pseudorapidity {rj). As noted in Section II. C, such 
a gap suppresses short-range correlations such as 
Bose-Einstein interference and Coulomb final-state 
interactions. To estimate the non-flow contributions 
to our estimate of 7j2{TPC} due to these short-range 
correlations, we have studied variations in the re- 
sulting tj 2 {TPC} induced by varying the event ver- 
tex selection along the beam direction, by varying 
the DCA cut value, and by varying the size of the 
pseudorapidity gap between the sub-events in the 77 
sub-event method. Tests of this type suggest that 
the systematic error on 7j 2 is on the order of 10%. 
However, based on an alternative approach [23| to 
fitting and interpreting the measured data at 62.4 
GcV and 200 GcV, it is argued that the systematic 
uncertainty could be larger. 



IV. SUMMARY 

In this study of Cu+Cu collisions at ^/snn = 22.4 
GeV, we present results at midrapidity and at for- 
ward rapidity for directed flow v\ and elliptic flow 
as a function of pscudorapidity, and for elliptic flow 
V2 as a function of transverse momentum at midra- 
pidity for charged hadrons. For our directed flow 
measurement, non-flow correlations are expected to 
be strongly suppressed by the large pseudorapid- 
ity gap between the detector used for event plane 
determination (BBC) and the main tracking detec- 
tor (TPC) used in our directed flow measurement. 
Our estimate of the systematic error on the di- 
rected flow 7)1 {BBC} is not more than 10% to 15% 
for the TPC and FTPC regions, respectively Our 
findings for Cu+Cu and Au+Au at 22.4 GeV ex- 



tend observations previously made at 62.4 GeV and 
200 GeV that directed flow is independent of sys- 
tem size in this energy region. Our findings also 
demonstrate that vi(i]/yb eam ) remains independent 
of beam energy down to 22.4 GeV, a scaling be- 
havior that has already been established at 62.4 
GeV and 200 GeV. We find that directed flow vio- 
lates the "entropy-driven" multiplicity scaling which 
dominates all other soft observablcs. An important 
feature associated with generation of v\ in the of 
the collision process is that different rapidity losses 
need to occur for particles at different distances from 
the center of the participant zone, which is beam- 
energy dependent. Measurements of the elliptic flow 
for 22.4 GeV Cu+Cu collisions are also presented. 
We compare pr-integrated V2 with measurements at 
higher energies. The pT-dcpcndence of the measured 
v 2 at 22.4 GeV is similar to that at 62.4 and 200 
GeV. UrQMD and AMPT models (the latter both 
with and without string melting) do not agree with 
the present measurements for both first (v±) and sec- 
ond (i>2) coefficients. 
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